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Abstract. First, we derive explicit computable expressions of structured backward errors of approx- 
imate eigenelements of structured matrix polynomials including symmetric, skew-symmetric, Her- 
mitian, skew-Hermitian, even and odd polynomials. We also determine minimal structured pertur- 
bations for which approximate eigenelements are exact eigenelements of the perturbed polynomials. 
Next, we analyze the effect of structure preserving linearizations of structured matrix polynomials 
on the structured backward errors of approximate eigenelements. We identify structure preserving 
linearizations which have almost no adverse effect on the structured backward errors of approximate 
eigenelements of the polynomials. Finally, we analyze structured pseudospectra of a structured matrix 
polynomial and establish a partial equality between unstructured and structured pseudospectra. 
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1 Introduction 

Consider a matrix polynomial P(z) := J2jLo zJ A? of degree m, where Aj S C nXTl and A m ^ 0. 
We assume that P is regular, that is, det(P(z)) ^ for some zeC. We say that A S C is 
an eigenvalue of P if det(P(A)) = 0. A nonzero vector x G C" (resp., y S C") that satisfies 
P(X)x = (resp., y H P(X) = 0) is called a right (resp., left) eigenvector of P corresponding to 
the eigenvalue A. The standard approach to computing eigenelements of P is to convert P into 
an equivalent linear polynomial L, called a linearization of P, and employ a numerically back- 
ward stable algorithm to compute the eigenelements of L, where L(z) := zX+Y, X 6 c mnxm " 
and Y € c«mx««\ j g we jj ]j nown that a matrix polynomial admits several linearizations. In 
fact, it is shown in [20|,I18| that potential linearizations of a matrix polynomial form a vector 
space. Thus choosing an optimal (in some sense) linearization of P is an important first step 
towards computing eigenelements of P. In general, a linearization of P can have an adverse 
effect on the conditioning of the eigenvalues of P (see, [IS])- Hence by analyzing the condition 
numbers of eigenvalues of linearizations, potential linearizations of P have been identified in 
[T5] whose eigenvalues are almost as sensitive to perturbations as that of P. Further, it is 
shown in [11) that these linearizations are consistent with the backward errors of approximate 
eigenelements in the sense that they nearly minimize the backward errors. 

Polynomial eigenvalue problems that occur in many applications possess some distinctive 
structures (e.g., Hermitian, even, odd and palindromic) which in turn induce certain spectral 
symmetries on the eigenvalues of the matrix polynomials (see, [2lJ H3 [24j [17l [16] and the ref- 
erences therein). With a view to preserving spectral symmetry in the computed eigenvalues 
(and possibly improved accuracy), there has been a lot of interests in developing structured 
preserving algorithms (see, [T21 [231 ESI US] and the references therein). Since linearization is 
the standard way to solve a polynomial eigenvalue problem, for a structured matrix polyno- 
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mial it is therefore necessary to choose a structured linearization and then solve the linear 
problem by a backward stable structure preserving algorithm. For the accuracy assessment 
of computed solution, it is therefore important to understand the sensitivity of eigenvalues of 
a structured matrix polynomial with respect to structure preserving perturbations. Also it is 
equally important to know the structured backward errors of approximate eigenelements of a 
structured matrix polynomial. Moreover, for a variety of structured polynomials such as sym- 
metric, skew-symmetric, Hermitian, skew-Hermitian, even, odd and palindromic polynomials, 
there are infinitely many structured linearizations, see |12[I21|. This poses a genuine problem 
of choosing one linearization over the other. For computational purposes, it is highly desirable 
to know how different structured linearizations affect the accuracy of computed eigenelements. 
Thus the selection of an optimal or a near optimal structured linearization is an important 
step in the solution process of a structured polynomial eigenvalue problem. The sensitivity 
analysis of eigenvalues of structured matrix polynomials with respect to structure preserving 
perturbation has been investigated in @J. It also provides a recipe for choosing structured 
linearizations whose eigenvalues are almost as sensitive to structure preserving perturbations 
as that of the structured matrix polynomials. 

To complete the investigation, in this paper we analyze structured backward errors of 
approximate eigenelements of symmetric, skew-symmetric, Hermitian, skew-Hermitian, T- 
even, T-odd, H-even and ii-odd polynomials. These structures are defined in Table [T] The 
main contribution of this paper is as follows. 

First, we derive explicit computable expressions for the structured backward errors of 
approximate eigenelements of structured matrix polynomials. We also construct a minimal 
structured perturbation so that an approximate eigenelement is the exact eigenelement of the 
structured perturbed polynomial. These results generalize similar results in [3] obtained for 
structured matrix pencils. 

Second, we consider structured linearizations that preserve spectral symmetry of a struc- 
tured matrix polynomial and compare the structured backward errors of approximate eigenele- 
ments with that of the structured polynomial. For example, a T-even matrix polynomial 
admits T-even as well as T-odd linearizations both of which preserve the spectral symmetry 
of the T-even polynomial. Based on these results we identify structured linearizations which 
are optimal in the sense that the structured backward errors of approximate eigenelements of 
the linearizations are bounded above and below by a small constant multiple of that of the 
structured polynomials. We show that these linearizations are consistent with the choice of 
linearizations discussed in [4] by analyzing structured condition numbers of eigenvalues. 

Third, we show that the effect of structure preserving linearization on the structured 
backward errors of approximate eigenelements is almost harmless for a wide class of structured 
linearizations. We show that bad effect, if any, of a structure preserving linearization can be 
neutralized by considering a complementary structured linearization. For example, when P 
is a T-even polynomial, we show that any T-even linearization is optimal for eigenvalues A 
of P such that |A| < 1, and any T-odd linearization is optimal for eigenvalues A such that 
|A| > 1. In such a case, we show that the backward error of an approximate eigenelement of 
the linearization differ from that of P by no more than a factor of 2. We show that similar 
results hold for other structured polynomials as well. In contrast, it is shown in [4] that 
the condition numbers of eigenvalues of these optimal linearizations differ from that of the 
polynomial by a factor of a constant whose size could of the order of the degree of the matrix 
polynomial. 

Finally, we analyze structured pseudospectra of structured matrix polynomials and estab- 
lish a partial equality between structured and unstructured pseudospectra. Similar study for 
palindromic matrix polynomials has been carried out in [2] , see also p] [8] . 

The rest of the paper is organized as follows. In section [2j we review structured poly- 
nomials and their spectral symmetries. In section [3j we analyze structured backward errors 
of approximate eigenpairs of structured polynomials. In section 31 we analyze the effect of 
structure preserving linearizations on the backward errors of approximate eigenelements of 
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structure polynomials and provide a recipe for choosing optimal linearizations. Finally, in 
section [5J we consider structured pseudospectra of structured matrix polynomials. 

2 Structured matrix polynomials 

We consider matrix polynomial of degree m of the form P(z) := YllLo ^Aj'' where Aj £ C nxn 
and A m ^ 0. Let P m (C" x ") denote the vector space of matrix polynomials of degree at most 
m. The spectrum of a regular polynomial P G P, m (C™ xrl ), denoted by <r(P), is given by 
<r(P) := {z E C : det(P(z)) = 0}. Strictly, speaking u(P) consists of finite eigenvalues of P. 
If the leading coefficient of P is singular then P has an infinite eigenvalue. In this paper, we 
consider only finite eigenvalues of matrix polynomials. An infinite eigenvalue of P, if any, can 
easily be analyzed by considering the reverse polynomial of P (see [6]). We say that (X,x,y) 
is an eigentriple of P if A is an eigenvalue of P and, x and y are the corresponding nonzero 
right and left eigenvectors, that is, P(X)x = and y H P(X) = 0. 

We denote the transpose and conjugate transpose of a matrix A by A T and A H , respec- 
tively. Define the map P m (C" xn ) -> P m (C" x "), P h-> P* given by P*(z) := £™ z^*, 

where A* = A T or A* = A H . The map P h> P* can be used to define interesting structured 
matrix polynomials such as symmetric, skew-symmetric, Hcrmitian, skew-Hermitian, *-even 
and *-odd matrix polynomials. These structures are defined in Table [TJ The table also shows 
the eigentriples as well as the spectral symmetries of the eigenvalues, see also [21]. We denote 
the set of structured polynomials having one of the structures given in Table [1] by §. By 
writing a pair (A, /x) in the third column of Table [T] we mean that if A is an eigenvalue of P 
then so is fj,. Notice that the eigenvalues of Hermitian and skew-Hermitian polynomials have 
the same spectral symmetry. Similarly, the eigenvalues of *-even and *-odd polynomials have 
the same spectral symmetry, where * G {T, H}. 



§ 


Condition 


spectral symmetry 


eigentriple 


symmetric 


P T (z) = P(z),VzeC 


A 


(A, x,x) 


skew-symmetric 


Y T {z) = -P(z),Vz e C 






T-even 


P T (z) = P(-z),Vz€ C 


(A, -A) 


(A,a;,f), (~X,y,x) 


T-odd 


P T (z) = -p(_ z ),Vze C 






Hermitian 


P H (z) = P(z),Vz G C 


(A, A) 


(X,x,y), (A, y, x) 


skew-hermitian 


P H (z) = -P(z),Vz G C 






-ff-even 


P H (z)=P(-z),VzeC 


(A, -A) 


(A, £,?/), (-X,y,x) 


iJ-odd 


P H (z) = -P(-z),Vz G C 







Table 1: Spectral symmetries of structured polynomials. 



Let P € § be regular. With a view to obtaining structured backward error of (A, x) G 
C x C™ with x H x — 1 as an approximate eigenpair of P, we now show that there always exists a 
polynomial AP G S such that (A, x) is a right eigenpair of P + AP, that is, (P(A) + AP(A))x = 
0. Recall that § denotes the set of structured polynomials having one of the structures given in 
TableQ] In short, we write S G {sym, skew-sym, Herm, skew-Herm, T-even, T-odd, iJ-even, ii-odd}. 

Theorem 2.1 LetE>& {sym, skew-sym, Herm, skew-Herm, T-even, T-odd, H -even, H -odd} and 

P G S be given by P(z) = Y^=o z3A j- Let ( X > x ) ^ C x C n be such that x H x = 1. Set 
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-P(A)a;, A m := [1, A, ... , A m ] T and P x := I - xx 11 . Define 



AA 3 



AA 3 



xx T AjXX H + ||A AJ ||a [xr T + rx H - 2(r T x)xx H ], if Aj = Af , 



'J^ I 2 \ xrT ~ rxH ^ if A i ~ A j> 

-xx H AjXX H + p^-jp [X>xr H P x + ~X?P x rx H ], if Aj = Af , 

-xx H AjXX H - ,, A \ 2 [X j xr H P x - >JP x rx H ], if Aj = -Af, 



and consider the polynomial AP(z) := Y^JLo z J AAj. Then P(A)x + AP(A)x = and AP e S. 
Proof: The proof is computational and is easy to check. ■ 

3 Structured backward errors 

Backward errors of approximate eigenelements of regular matrix polynomials have been sys- 
tematically analyzed and computable expressions for the backward errors have been derived 
by Tisseur in [26] . For our purpose, we require a different norm setup for matrix poly- 
nomials. We equip V m (C nxn ) with a norm so that the resulting normed linear space can 
be used for perturbation analysis of matrix polynomials. Let P S P m (C™ xn ) be given by 

p ( z ) == Ejlo A i zj - We define 

1/2 

.1 w | W A o\\m 

where ||A||m denotes the Frobenius norm when M = F and the spectral norm when M = 
2. Accordingly, we say that \\-\\f is the Frobenius norm and |||-|||2 is the spectral norm on 
P m (C raXTl ). See [6l[5] for more on norms of matrix polynomials. 

Let (A, x) S C x C" be such that x H x = 1 and P e P m (C" xn ) be regular. We denote the 
backward error of (A, x) as an approximate eigenelement of P by t]m(\ x, P) given by 

r) M (\,x,P):= inf {||| AP||| M : P(A)z + AP(X)x = 0}. 

APeP m (t> x ' 1 ) 

Setting r := — P(A)x and A m := [1, A, ... , A m ] T , it is easily seen that 

A-^ TX^ 

for M — F as well as M — 2. Indeed, defining A Aj := „ — j = : m, and considering 

the polynomial AP(z) := jy™ =0 z-* AAj, we have |||AP|||m = H^lb/H^IUHAml^ and P(A)a; + 
AP(A)x = 0. Consequently, for simplicity of notation, we denote 7]m{X, x, P) by ry(A,a;,P). 

Now suppose that P e §. Then treating (A,x) as an approximate eigenelement of P, we 
define the structured backward error of (A, x) by 

V S M (\x,P) := A inf s {l|AP||| M :P(A)a ; + AP(A)x = Q}. 

In view of Theorem |2~TI it follows that r?(A,x,P) < nlj(X,x,P) < oo. Structured backward 
errors of approximate eigenelements of structured matrix pencils have been systematically 
analyzed and computable expressions of the structured backward errors have been derived 
in [3]. In this section we generalize these results to the case of structured matrix polynomials. 
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As we shall see, determining 77^ (A, x, P) is much more difficult than determining rj%(X, x, P) 
and requires solution of norm preserving dilation problem for matrices. The Davis-Kahan- 
Weinberger solutions of norm preserving dilation problem given below will play an important 
role in the subsequent development. Let A, B, C and D be matrices of appropriate sizes. 
Then the following result holds. 



Theorem 3.1 (Davis-Kahan- Weinberger, [10J ) Let A,B,C satisfy 



[i and 



A C 
B D 



|| [A C] || = /i. Then there exists D such that 
this property are exactly those of the form 

-KA H L + nil - KK H ) l ' 2 Z{I - L H L) 1 ' 2 , 



pi. Indeed, those D which have 



D 

where K H := ([i 2 I—A H A)~ 1 / 2 B H , L := {pi 2 l— AA H )~ 1 / 2 C and Z is an arbitrary contraction, 
that is, \\Z\\ 2 < 1. ■ 



For a more general version of the above result, see [ID] . 



3.1 Symmetric and skew-symmetric polynomials 

We now derive structured backward error of (A, x) € C x C™ as an approximate eigenpair 
of symmetric and skew-symmetric matrix polynomials. We also derive minimal structured 
perturbations so that (A, x) is an exact eigenpair of the perturbed polynomials. First, we 
consider symmetric matrix polynomials. Note that a matrix polynomial P € P m (C" xn ) is 
symmetric if and only if all the coefficient matrices of P are symmetric. For a symmetric 
matrix polynomial, we have the following result. 

Theorem 3.2 Let § denote the set of symmetric matrix polynomials in P m (C nx ") and let 
Pe§. Let (X,x) eCxC" be such that x H x = 1. Set r := -P(X)x, P x := I — 



xx H and 



A m := [1, A, ... , X m ] T . Then we have 



r,%(X, x, P) = |xTr|2 - < V2n{X, x, P) and rg (A, x, P) = ^(A, x, P). 

Am 2 



^ J ?[xr T + rx H — (r T x)xx H ], j — : to, and consider the polynomial AP(z) 



Set AAj : 

Y^jLo AAj. Then AP is a unique polynomial such that AP G S, AP(A)i + P(A).r = and 
|||AP|||i? = r^r(X, x, P). Further, define 



[xr + rx — (r x)xx 



I -Am || 2 



Xi x T r PTrr T P^ 



lAmlll 



|a; T r|^ 



ana consi 



IIAPH 



der the polynomial AP(z) := 52™ =0 AAj . Then AP e S, AP(A)i + P(A)a: = and 
vt(X,x,P). 



Proof: In view of Theorem EU let AP € S given by AP(z) := ^JL AAj-z^ be such that 
P(X)x + AP(A)x = 0. Let Qi e C™*^-^ be such that the matrix Q = [x Q{\ is unitary. 



Then AAj := Q T AA 3 Q = 
we have 



'■]} 



3 

A, 



, where Xj — Xj is of size n — 1. Since QQ T = I, 



Q(AP{X))Q H x = r^> (AP{X))Q H x = Q 1 r = 



T 

x r 
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T 

x r 



( Y\ _ X J a ■ 

As Q H x = ei, the first column of the identity matrix, we have I -,—jm \,- J3 \ — \ 

_ V l^j=a A a j J V 

Hence the minimum norm solutions are aj — y^. |p an< ^ a n = \\\ X jp 1 j = : m. Conse- 
quently, we have 



(||A„il 5 



(2) 



of A 

\AAj\\ 2 F = \\AAj\\ 2 F = |a^| 2 + 2||a J || 2 \ Since C^Qf = I -xx T , we have 



This shows that the Frobenius norm of AAj is minimized when Xj = 0. Hence we have 



Now from we have 



-P , 2||(/-x^)r||| _ y/2\\ 



r\\l 



\h-m\\ 2 2 || A-m.|| 2 ||A m ||2 



\ix T r 




l|A m ||l 

which gives the desired polynomial AP for the Frobenius norm. 

For the spectral norm, we employ dilation result in Theorem 13.11 to the matrix in Q . 

Indeed, for ^ := 1 4^pK by Theorem O we have 



AJ x?r Qfr(Q{r) T ._ 
J IIAmllKIHIl-l^rP)^ U ' W ' 

which gives ?7§(A, x, P) = Jj ? = rj(X, x,P). Putting Xj in ^ and after simplification we 

|| A m ||2 

have 



AAj = [xr T + rx H - (r T x)xx H ] - 



\Am\\r v ' J ||A m ||l(||r||l-|^r|2) 

which gives the desired polynomial AP for the spectral norm. ■ 

Remark 3.3 If \x T r\ = ||r|| 2 , then ||Q^r|| 2 = 0. Hence considering Xj = 0, j = : to, in 
the above proof we obtain the desired results for the spectral norm. Note that in such a case 
we have rf F 

Observe that if Y is symmetric and Yx — then Y — P x ZP X for some symmetric matrix 
Z. Consequently, from the proof Theorem 13.21 we have QjXjQ^ = Pf?ZjP x ,j — : m, for 
some symmetric matrices Zj. Hence we have following. 

Corollary 3.4 Let P be a symmetric matrix polynomial. For (A, x) G C x C" with x H x = 1, 
set r := — P(A)a;. Then there is a symmetric matrix polynomial Q such that P(A)a;+Q(A)iz; = 
if and only ifQ(z) = AP(z) + P^R(z)P x for some symmetric polynomial R, where AP is the 
symmetric polynomial given by AP(z) := 53j=q 2,3 ^A? an d 



AAj := AJ ||2 [xr T + rx H - (r T x)xx H ], j = 



m 1 1 2 



Next, we consider skew-symmetric matrix polynomials. Note that a matrix polynomial 
P € P m (C™ x ") is skew-symmetric if and only if all the coefficient matrices of P are skew- 
symmetric. For skew-symmetric matrix polynomials we have the following result. 
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Theorem 3.5 Let § denote the set of skew- symmetric matrix polynomials in P m (C nx ") and 
let Pe§. For (A, x) G C x C" with x H x = 1, set r := -P(A)x. Then we have 

rf F {\x,V) = V2 r,(X,x,P), r|(A,z,P) = ??(A,z,P). 

For the skew- symmetric polynomial AP given in Theorem \2.1\ we have P(A)a? + AP(A)x 
0. lAPlli? =7?|(A,x,P) and |||AP|||2 = rjf(A, a;, P). 



Proof: The proof is the same as that of Theorem 13.21 except that AAj is skew-symmetric for 
j = : m. This gives 



-( 



l|A m || 



A^=| Jj QTr ^ I • (3) 



Setting = 0, we obtain the results for the Frobenius norm. 

Setting Hj := ^^JJ ^l 2 and invoking Theorem 13. li it is easily seen that the spectral norm 

of AAj in ([3]) is minimized when Xj = 0. Hence the desired results follow for the spectral 
norm. ■ 

Note that if Y is a skew-symmetric matrix and Yx = then Y = P x ZP X for some 
skew-symmetric matrix Z. Hence we have the following result. 

Corollary 3.6 Let P G P m (C nx ") be a skew- symmetric matrix polynomial. For (A, x) G C x 
C n with x H x = 1, set r := — P(A)x. Then there is a skew- symmetric matrix polynomial Q such 
that Y(X)x + Q(\)x = if and only if Q(z) = AP(z) + P^K(z)P x for some skew- symmetric 
polynomial R, where AP is the skew- symmetric polynomial given by AP(z) := ^™L z l AA, 

and 

X> 
I A, 



AAj = - n f 2 [xr T - rar 8 ], j = : m. 



Mn 1 1 2 



3.2 T-even and T-odd matrix polynomials 

For backward perturbation analysis of T-even and T-odd polynomials, we need the even index 
projection n e : C m+1 — > C m+1 given by 

., , T _ J [x , 0, x 2 , 0, . . . , x m - 2 , 0, x m ] T , if m is even, 

n eU a;o, a*, ... , sm_i, x m j j .- | ^ q, ar 2 , 0, . . . , 0, x m _i, 0] T , if m is odd. 

Note that "0" is considered as even number. Observe that I — H e is the odd index projection. 

Recall that a matrix polynomial P G P m (C™ xn ) given by P(z) := J2jLo ^i z ^ xs T-even if 
and only if Aj is symmetric when j is even (including j = 0) and Aj is skew-symmetric when 
j is odd. We have the following result for T-even matrix polynomials. 

Theorem 3.7 Let § denote the set of T-even matrix polynomials in P m (C" x "). Let P G § 
and (X,x) G C x C" be such that x H x = 1. Set r := -P(X)x, P x := L 
[1, A, ... , A m ] T . Then we have 



\x T r\ 2 , J|r||l - |^ r |a \ x T r? ||r||I - |* T rf 



ri%(X,x,P) = x l ^ - ' +2 " m . ' 1 , 4(X,x,P)- v , 



|n e (A ro )||i ||A m ||^ ' ' ' > W ||n e (A m )||^ ||A 



mil 2 



In particular, if m is odd and |A| = 1 i/ien we focwe rj F (X, x.P) = a/2 7y(A,x,P) and 
r ? f(A,x,P)=r 7 (A, a; ,P). 
Tor j = : m, define 



^{x T r)xx H + — - — [xr T Pa: + Tjrx ff ], i/j is even, 



|rie(A m )||l^ ' ||A m ||| 



|A m ||| 



P x rx H — xr T P x ], if j is odd . 
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Then AP(z) := Ei=o z " ' * s a unique T-even polynomial in § such that P(A)x + AP(A)i = 
and |||AP|||i? = ?7p(A,x,P). Further, for j = : m, defining 

)J Wr P^rr T P x 



AAi 



E 



* ||n e (A m )||l (||r||! - |^rf )< 



i/j is even, 



^ , if j is odd, 

we obtain a T-even polynomial AP(z) := Ej— p z J AA,- in § such that P(A).t + AP(A).t = 
and }APfo=r${\,x,P). 

Proof: In view of Theorem 123] let AP e § be such that P(X)x + AP(A)i = 0. Assum- 
ing that AP is given by AP(z) := Ei=o £^AjZ 3 , and arguing similarly as in the proofs of 



Theorems 13.21 and 13.51 we have A A 



a rj aj 



J i I . Xj = Xj when j is even, and 

3 ^3 / 



AA 3 = 



bj 

~bj Yj 



Yj = —Yj when j is odd. Consequently we have 



Ej x '<>jj 
Ej-even X ° a 3 - Ej-odd Xlb 3 



T 

x r 



Hence the smallest norm solutions are a 
Therefore, we have 



33 - ||n e (A m )|| 2 xTr ' a 3 - \\A m \\'iQl r > b 3 ~ Pvfllf ^l 7 "* 



\i T 

\\n c (A m )\\ 2 2 x r MIA 



AAj = Q T AAjQ = { 



X,; 



if j is even 



(4) 



I V 



if j is odd. 



Setting Xj = = Yj and using the fact that QiQ^ = I — a;a; T , we obtain the desired unique 
T-even polynomial AP(z) := Ej=o z ^Ej such that 



||AP||| F =r ? |(A,a:,P) = 



I 21 ' 7,7 '! 2 o H r ll2 — | a;Tr | 2 



|n e (A m )| 



I A TO || 2 



When to is odd and |A| = 1, it is easily seen that ||n e (A m )||| = |||A m |||. Hence we have 
7/|(A,a;,P) = y/2 v (\,x,P). 



/|\j|2 \ x T r \2 lAJ'l 2 dlrll — \x T r\ 2 ) 

For the spectral norm, setting fj,j '■= y Wj ^ )|| 4 + — — a !|| 4 when j is even, and 

■ — ■ — . 2 J — when j is odd, and applying Theorem 13.11 to the matrices in HJ), 

||A m || 2 

AJ Wr Ql r{Ql r) T 



fij := 
we have 



X > ||n e (A m )||l (Hrlll-l^rl 2 ) 



and Yj = 0. 



Consequently, we have r?|(A, x, P) 



|a; T r| 2 



|n e (A m )||^ \\A m \\i 



From we have 



AAt 



SrJT^Li + liTTia ^ rTp - + P - Trxff 1 + Qi x iQ?y if J is even 
||li e (A m J|| 2 ||A|| 2 



A' 



if j is odd. 
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Substituting Xj and Yj in AAj we obtain the desired T-even matrix polynomial AP for the 
spectral norm. ■ 

Remark 3.8 If\x T r\ = ||r||2 then ||<3^V||2 = 0. Hence considering Xj = = Yj in the above 
proof, we obtain the desired result for the spectral norm. Note that in such a case we have 
if F (X,x,P) - V2ri(\x,P) = V2rj(X,x,P). 

Recall that when A is symmetric (resp., skew-symmetric) and Ax = then A = P x ZP X for 
some symmetric (resp., skew-symmetric) matrix Z. Consequently, from the proof Theorem [3~71 
it follows that AA, := Ej + P x ZjP x , where Zj = Zj when j is even, and Zj = —Zj when j 
is odd. Hence we have the following result. 

Corollary 3.9 Let P be a T-even matrix polynomial in P m (C™ xn ). Let (X,x) G C x C" be 
such that x H x = 1. Then there is a T-even matrix polynomial Q such that P(A)x + Q(A)x = 
if and only i/QO) = AP{z) + P x r R(z)P x for some T-even matrix polynomial R G P TO (C" X ™), 
where AP(z) :— 'Y^^EjZ- 1 and Ej 's are given in Theorem \3.7\ 

Next, we consider backward error of T-odd polynomials. Observe that a matrix polynomial 
P G P m (C nx ") given by P(z) := J^jLo ^o z ^ ^ s ^-°dd if and only if Aj is skew-symmetric when 
j is even (including j — 0) and Aj is symmetric when j is odd. 

Theorem 3.10 Let § denote the set of T-odd matrix polynomials in P TO (C" X ™). Let P G §> 
and (X,x) G C x C" be such that x H x = 1. Set r := — P(A)x, P x := I — xx H and A m := 
[1,A,...,A 




|(j-n c )(A m )||^ + 1 — \\A^\p 2 — ' y A ^ u > 
\^r?(A,^P), " ifX = 0, 



\x T r\ 2 , \\r\\l-\x T r\ 2 ., , , „ 

||(/-n e )(A m )||^ + ||A m ||3 ' 
T](X,x,P), ifX = 0. 



In particular, ifm is odd and \X\ — 1 we have ij^(X, x, P) = y/2 r)(X, x, P) and ?7|(A, x, P) 
n(X,x,P). For j — : m, define 



[P x rx H — xr T P x ], if j is even 



I A II 2 



X^xx T rx H A J 
\\(I-U e )(A m )\\l + JX mll2 



2 + ttt — [xr T P x + F x rx H ], if j is odd. 



Then AP(z) := 53j=o z " ' Fj is a unique T-odd polynomial in S such that P(X)x + AP(A)x = 
and |||AP||| F = tj%(X,x,P). 

Further, for j = : m, define AAj :— Fj when j is even, and 



a a =F AJ x T rP^rr T P x 

J ' J ||(7-n e )(A m )|||(||r||I-|^rP) 

when j is odd. Then AP(z) :— X)J=o zJ AAj is a T-odd polynomial in § such that P(A)a; + 
AP{X)x = and \\AP\l 2 = ?jf (A, x, P). 

Proof: The desired results follow from the proof of Theorem 13. 71 by interchanging the role of 
AAj for even j and odd j. I 

We have the following results whose proof is immediate. 

Corollary 3.11 Let P be a T-odd matrix polynomial in P m (C" x "). Let (A, x) G C x C" be 
such that x H x = 1. Then there is a T-odd matrix polynomial Q such that P(X)x + Q(X)x = 
if and only if Q(z) = AP(z) + P x K(z)P x for some T-odd matrix polynomial R G P m (C nx "), 
where AP(z) := X)^=o Fj% 3 and Fj 's are given in Theorem \3.10\ 
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3.3 Hermitian and skew-Hermitian matrix polynomials 

We now consider structured backward errors of approximate eigenelements of Hermitian and 
skew-Hermitian matrix polynomials. We proceed as follows. Let S C P m (C" xrl ) and uj G C 
be such that |w| = 1. We set := {ujP : P G §}. Then for P € P m (C nx "), it is easily seen 
that 

tj|(A, x, P) = rfj? (A, x, wP) and r?f (A, x, P) = r%» (A, x, wP). (5) 

Note that a matrix polynomial P G P m (C" xn ) is Hermitian (resp., skew-Hermitian) if and 
only if all the coefficient matrices of P are Hermitian (resp., skew-Hermitian). Let Herm and 
skew-Herm, respectively, denote the set of Hermitian and skew-Hermitian matrix polynomials 
in P m (C" xn ). Then noting that a matrix X S C" xn is Hermitian if and only if iX is skew- 
Hermitian, it easily seen that the maps 



Herm — > skew-Herm, P i — > iP and skew-Herm — > Herm, Q i 



(6) 



are isometric isomorphisms. Thus, in view of ([5|) and ((6|), it follows that the structured 
backward error of (X,x) as an approximate eigenpair of a skew-Hermitian polynomial can 
be obtained from the structured backward error of (A, x) as an approximate eigenpair of 
a Hermitian matrix polynomial and vice-versa. We therefore analyze structured backward 
perturbation of Hermitian matrix polynomials. 

For i € C™, we denote by Re (a;) and lm(x), respectively, the real and the imaginary parts 
of x. Then we have x = Re(a;) 4- ilm(a;). We denote the real and imaginary part of a complex 
number z G C by re(z) and im(z), respectively. We denote the Moore-Penrose pseudo-inverse 
of A by and the canonical basis of C m+1 by ej,j = : m. 



Theorem 3.12 Let Herm denote the set of Hermitian matrix polynomials in P m (C nx "). Let 
P e Herm and (A, x) G C x C™ be such that x H x = 1. Set r := -P(A)x, P x := I 
A m := [1, A, ... , A m ] T . Then we have 



xx H and 



^ erm (A,x,P) 



Herm (A,x,P) 



-\x H r\ 2 



2(lkil 



< V2r)(\ ,x,P), 

f-\x H r\ 2 ) 



r)(\,x,P), 



i | r j [ | — | a: 77 r | 2 



z/Ae 
z/Ae 



i/AeK, 
if A g C \ 

\M, 



where r 



"Re(A m ) T " 


t 


re(x H r) 


lm(A m ) T _ 




\m(x H r) 



. For the Frobenius norm, define 



X' 



AA, 



l|A„..._ 
ejrxx + 



;(xr H + rx H — (r x)xx H ), 



when A G 



|A m || 2 



[\3P x rx H + X>xr H P x ], when\eC\\ 



Then AP(z) := X)J=o z ^ * s a unique Hermitian polynomial in Herm such that P(A)a; 
AP(A)a; = and |||AP||| F = ^ erm (A, x, P). 
for the spectral norm, define 



A-' 



A An 



|A m ||| 



(rx ff + .ir fl — (r x)xx H ) — 



X> x H rP r rr H P r 



IA II 2 



<n|l2 



[A-J'P^ra^ + Xxr H P x ] - / 



2 _ 1^2)' 

eTr P x rr H P T 



(l|r|| a 



r s — ar J r 



ff».|2 ' 



w/iew A G K, 
when A G C \ . 



TTiera AP(z) := X)J=o zJ AA,- * s a Hermitian polynomial in Herm sucft i/iat P(A)ar+AP(A)a; 
and |||AP||| 2 = r/^ erm (A,x,P). 
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Proof: Again, in view of Theorem l2.11 let AP(z) = Y^j=o z ~* be a Hermitian polynomial 
such that AP(A)i + P(X)x = 0. Choosing a unitary matrix Q :— [x, Qx], we have 

/ y"?L A?a- \ / x H r \ 
Now AP(A)i + P(A)a; = u " = n H ■ The minimum norm solution of 

Eilo AJfl i = Q? r is s iven b y a J = jx^Qi r - 

Now suppose that A £ M. Then the minimum norm solution of Y^JLo ^ a ii = ^i^ 7 " is given 
by Qjj = || A AJ |p x H r 6 M. Hence for A 6 K, we have 



~W X r V ||A„, ill ^1 r J 



u /2||r|| 2 — Ir-^a;! 2 

For the Frobenius norm, setting Xj = we obtain r/^ erm (A, x, P) = v | ^ ^ — — and the 



For the spectral norm, setting /ij := ^'J 2 and applying Theorem 13. II to (|7|), we obtain 



desired Hermitian polynomial AP 

„ — P 

|| TV 

X j x H r{Q^r)(Q^r) H 
j ||A m |||(||r|||-|^r|2)- 

This gives rj Herm (\,x,P) = ||Pj| 2 = ^(A, x,P). Now substituting Xj in and simplifying 
the expression, we obtain the desired Hermitian polynomial AP. 

Next, suppose that A e C \ K. Then the minimum norm solution of X)J=o ^ a jj = x r * s 
obtained by solving 



Therefore we have ajj — ejr. Hence for A G C \ E, we have 



Q H AA j Q=[ A / J ' „ V||A ™"= - 1 , j = 0:m. (8) 



Thus, for the Frobenius norm, setting Xj = we obtain 



I A II 2 



and the desired Hermitian polynomial AP. 

For the spectral norm, setting fj,j := ^J\ejr\ 2 + ^ ^iia 2 ^ an< ^ a PPl vm g Theo- 
rem [Jl] to the matrix in |8]), we have 



ejf (Qf r)(Qf r) g 
\x H r\ 

This gives 

??2 Herm (A,x,P) = W||f 



|A m ||| 
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Now substituting Xj in ([5]) and simplifying the expression, we have 



1 rT -„ H m T-i i eJrP x rr H P x 



AAj = efrxx H + , lt ..J \iP x rx H + X J xr H P x 



A m ||I L *> \\ r \\l 



Hence the results follow. 



Remark 3.13 If \x H r\ = \\r\\2 then \\Qi -r || 2 = 0. Hence considering Xj = 0, j = : m, we 
obtain the desired results for the spectral norm. 

Let x £ C n be such that x H x = 1. If v4 E C nx ™ is Hermitian and Ax = then it is easily 
seen that A = (I — xx H )Z{I — xx H ) for some Hermitian matrix Z. Consequently, in view 
of Theorem 13.121 we have an analogue of the result in Corollary 13.41 for Hermitian matrix 
polynomials. 

Note that, in view of (© and ([6|), structured backward error of (X,x) as an approximate 
eigenpair of a skew-Hermitian matrix polynomial follows from Theorem 13.121 Indeed, let 
Q £ skew-Herm C P m (C™ x ") be a skew-Hermitian matrix polynomial. Then P := iQ £ 
Herm C P m (C" x "). Hence by © and ©, we have ^ ew_Herm (A, x, Q) = ?7^ erm (A, x, P). Now, 
let AP be the matrix polynomial given in Theorem 13.121 such that P(X)x + AP(A)x = and 
|||AP|||m = ?7M rm (A, x,P). Then setting AQ := -iAP, we have AQ £ skew-Herm such that 
Q(X)x + AQ{X)x = and ||| AQ| M = ^ ew - Herm (A, x, Q). 



3.4 H-even and H-odd matrix polynomials 

We now derive structured backward errors of approximate eigenelements of H-even and H- 
odd matrix polynomials. Recall that a matrix polynomial P £ P m (C™ xn ) given by P(z) := 
y^-o AjZ 3 is H-even if and only if Aj is Hermitian when j is even (including j = 0) and 
Aj is skew-Hermitian when j is odd. Let H-even and iJ-odd, respectively, denote the set 
of iJ-even and H-odd matrix polynomials in P m (C nx "). Then, as in the case of Hermitian 
matrix polynomials in ([5]), it is easily seen that the map 



H-even — ► H-odd, P 1 — ► iP and _ff-odd — > H-even, Q 



iQ 



(9) 



are isometric isomorphisms. Consequently, we only need to prove the results either for H-even 
or for H-odd matrix polynomials. Recall that A^ is the Moore-Penrose pseudo-inverse of A 
and ej,j — : m, is the canonical basis of C m+1 . 



Theorem 3.14 Set S := H-even C 

x H x 



:™ xn ). Let P £ § and (X,x) eCxC" be such that 



1. Set r := —P(X)x, P x := I — xx H and A m := [1, A, ... , X m ] T . Then we have 



where f 



4(x,x,p) = 



n e Re(A m ) T 
n e lm(A m ) T 



-\x H r\ 2 



<V^v(\x,P), ifXei] 



2(\\r\\i-\x»rn 



ifXe C\iR, 



T)(\,X,P), 



if X £ i R, 



(J-n e )lm(A ro ) 
(I-n e )Re(A m ) T 



A 

Tit 



^fii!, ifxeC\u 



re(x H r) 
\m(x H r) 



For j — : m, set 



E, 



' <XiP x rx H 



|A m ||l 



X°xr H P x ] and Fa 



! [XiP x rx H - X j xr H P x ]. 



I A II 2 
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For the Frobenius norm, define AAa := -rri — + r x H — (r x)xx H ] when A £ iK, and 

' IIAmlll 

J ejrxx 11 + Ej, if j is even, 
1 ' \ ieJrxx H + Fj , if j is odd, 

when AeC\il, for j = : m. Then AP(z) := EJlo z J AAj is a unique H-even polynomial 
in § such that P(X)x + AP(A)i = and ||| AP||| F = rfp(\, x, P). 
For the spectral norm, define 

~>J H H H H] (-l) 3+1 \ 3 x H rP x rr H P x 
r '~jKM [Xr +rX X)XX J+ l|A ro ||I(||r||l-|^rP) 

when A € iR, and 



, (-iy +1 ejr P x rr H l x 
e,- rxx + Ej + , _ , «/ j ewen, 



2 

* (-l ) i+1 ejf P^ P, 
rllo — \x H r\- 



iejrxx 11 + Fj n„n2 3 i [2 ~ ' * s oc ^> 



w/ien A £ C\ iM. TTien AP(z) := Ej=o ^'^A? * s an H-even polynomial in § smc/i i/iai 
P(X)x + AP{X)x = and |||AP||| 2 = r?| (A, z, P). 

Proof: By Theorem 12.11 there exists an H-even matrix polynomial AP(z) = Ej=o z,7 ^A; 
such that AP(A) = r. Now choosing a unitary matrix Q := [x, Qi], we have AAj = 

a " f ) Q H 1 Xf = Xj if j is even, and AAj =q( ™ rj °f \ Q H , Yf = —Yj 
a j Aj J \ a j fj J 

if j is odd. Notice that ajj is real for all j. 

Then AP(A)x = r gives ( %evea ^ + < Ej-odd ) = f ^ \ The mini _ 

V Eleven AJa j-£j-odd A a J / I r /_ 
mum norm solution of ^j-even ~ Ej-odd = Q? r is g iven by a j — \\kL\\1 ^ t ^ 3 
is even, and aj = ~jx^-jiQi !r if j 1S odd. 

Now suppose that A £ iR. Then the minimum norm solution for ctjj is given by ajj = 
rpAp X r when j is even, and eijj = — ^ x '^ x H r when j is odd. Hence ajj £ M when j is 
even, and iajj £ i K when j is odd. Consequently, we have 

_ " (10) 

II Am III ^1 r / 

(11) 

AJ Qfr 




when j is odd. Setting Xj = = Fj in (JTDJ) and (33]), we obtain ?j|(A,a;,P) = ^^aJ^' - 
and the desired AAj. 

Next, suppose that A e C \ i K. Then Ej-even •^"' a ji + * Ej-odd A5a ±; = xHr gi yes 

Ej-even re ( AJ ) a M-Ej-odd im ( A ^%J \ ( re{x H r) 
Ej-even im ( Ai ) a w + Ej-odd re ( AJ ) a J"i / V Mx H r) 
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Hence we have 



an 



V 



n e Re(A m ) T -(/-n e )lm(A m ) T 
n e lm(A m ) T + (/-n e )Re(A m ) T 



re(x H r) 
\m(x H r) 



= r^ ajj = e j r. 



J 



Consequently, we have 

Q H AAjQ 

when j is even, and 



Q H AAjQ 



when j is odd. Now setting Xj = 

AAj, j = : m, and ri%(\, x, P) = 
Frobenius norm. 



e 3 r 



A3 



x, 



(12) 



zej r 



Aj , 



(13) 



iiA m in 

= 1} in (fT2|l and ([13]), we have the desired matrices 



\x H r\' 



|A m ||l 



This completes the proof for the 



For the spectral norm, consider fij := when A G jR. Then applying Theorem 13.1 

to the matrices in (fTQ|) and (fTTj) . we obtain 

\i x H r{Q^r){Q^r) H 



X, 



, Y = V x H r(Q?r)(Q»r) H 

||A ro ||I(||r||2-|^P) and ' ||A m ||HIIHIi-|^l 2 )' 



Ik || 2 



This gives fg(A,x,P) - 
desired matrices AA,-, j = : m. 

When A G C \ i R, considering fij : 
rem 13. II to the matrices in (TT21 and (p~5|) . we obtain 



Now substituting X, and Yj in (JTUJ) and (HH), we obtain the 



r||| — k H r| 2 ) 



A m ||| 



and applying Theo- 



ejr (Qfr)(gfr) ff 

X J = - ii i — 



Consequently, we have 



\x H r\< 



^(A,z,P) 



and Y} = — 



iejf (Qf r)(Qf r) g 
llrllo - Ix-^rl 2 



|a; H r| 2 



|A m ||| 



Substituting Xj and K,- in (fT2"| and (jT3j) , we obtain the desired matrices AAj, j = : m. ■ 

Let i £ C be such that x^ir = 1. If X € C nX ™ j s skew-Hermitian and Xa; = then it is 
easily seen that X = (I — xx H )Z{I — xx H ) for some skew-Hermitian matrix Z. Consequently, 
it follows that an analogue of the result in Corollarv l3. 1 ll holds for H-even matrix polynomials. 

Observe that, in view of ([5]) and (O, the structured backward error of (A, x) as an 
approximate eigenpair of an _ff-odd matrix polynomial follows from Theorem 13.141 In- 
deed, let Q be an ff-odd matrix polynomial in P m (C nx "). Set § e := H-even C P m (C" xn ) 
and § := ff-odd C P m (C nx "). Then P := iQ e S e . Hence by © and ©, we have 
T)^(X, x, Q) = 7?m(A, x, P). Now, let AP be the matrix polynomial given in Theorem l3 . 141 such 
that AP e S e , P(A)at + AP(A)x = and |||AP||| M = rf^{\, x, P). Then setting AQ := -iAP, 
we have AQ e § , Q(A)x + AQ{\)x = and |||AQ||| M = rfy(\,x,P). 
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3.5 Polynomials with coefficients in Lie and Jordan algebras 



We mention that the structured backward perturbation analysis of structured matrix polyno- 
mials discussed so far can easily be extended to more general structured matrix polynomials in 
which the coefficient matrices are elements of appropriate Jordan and/or Lie algebras. Indeed, 
let M be a unitary matrix such that M T = M or M T = —M. Consider the Jordan algebra 
J := {A G C nxn : M~ 1 A T M = A} and the Lie algebra L := {A G C" x " : M~ 1 A T M = -A} 
associated with the scalar product (x, y) i— > y T M x. Consider a polynomial P(z) := Sj=o z3 Aj ■ 
Then by imposing the condition that the polynomial MP given by MP(z) = jyjLo X-'MAj is 
either symmetric or skew-symmetric or T-even or T-odd, we obtain various structured matrix 
polynomials. Said differently, S C P m (C™ xn ) defines a class of structured matrix polynomials 
if MS G {sym, skew-sym, T-even, T-odd}. Hence if P G § then the results obtained in the 
previous section are easily extended to P by replacing Aj and r := —P(X)x by MAj and Mr, 
respectively. 

Similarly, when M is unitary and M — M H or M — — M H , we consider the Jordan 
algebra J := {A G C nxn : M~ 1 A H M = A} and the Lie algebra L := {A G C nxn : 
M~ 1 A H M = —A} associated with the scalar product (x, y) i— > y H Mx. Then a class of 
structured matrix polynomials § C P TO (C™ X ™) is obtained by imposing the condition that 
MS G {Herm, skew-Herm, H-even, iTodd}. Hence the results obtained in the previous section 
are easily extended to P € S by replacing Aj and r 

In particular, when M := J, where J := 



-I 



= — P(A)x by MAj and Mr, respectively. 
€ C 2 " x2n , the Jordan algebra J consists 
consists of Hamiltonian matrices. So, for 



of skew-Hamiltonian matrices and the Lie algebra 
example, considering the polynomial P(z) := X)J=o zJ ^j, where AjS are Hamiltonian when j 
is even and skew-Hamiltonian when j is odd, we see that the polynomial JP(z) = X)j=o JAj 
is H-even. Hence extending the results obtained for iJ-even polynomial to the case of P, we 
have the following. 



Theorem 3.15 Let § denote set of polynomials of the form P(z) — 'Y^T = QZ 3 Aj where Aj 
is Hamiltonian when j is even, and Aj is skew-Hamiltonian when j is odd. Let P G § 
and (\,x) G C x C" be such that x H x = 1. Set r := -P(A)x, P x := I - 
[1, A, ... , A m ] T . Then we have 



xx H and A r 



V%(X,x,P) 



% S (A,x,P) = 



y/2\\r\\l-\xHjr\* 
l|A m || 2 



< V2r,(\,x,P), z/AG 



,2 , 2(||r||^-|x^./rP) 
Kilo i — !T2 : 



2 n ||A m ||| ' if X € C\i] 

n(X,x,P), ifXeiR, 

||r||^-|a;»JrP 



+ 



-, if X e C \ iM, 



whe 



n e (Re(A m ) T )-(/-n e )(lm(A m ) T ) 
n e (lm(A m ) T ) + (/-n e )(Re(A m ) T ) 



1 T re{x H Jr) 
\m(x H Jr) 



4 Effect of structured linearization on backward error 

As we have mentioned before, linearization is the standard approach to solving a polynomial 
eigenvalue value problem. It is well known that important classes of structured matrix poly- 
nomials admit structured linearizations jT5J Q21 US HO] • However, the process of linearizing a 
matrix polynomial (structure preserving or not) has its side effect too. It increases the sensi- 
tivity of eigenvalues of the matrix polynomial (see, [T31[TJ[3]). Therefore, it is important to 
identify linearizations whose eigenelements are almost as sensitive to perturbations as those of 
the matrix polynomial. Obviously, condition numbers of eigenvalues and backward errors of 
approximate eigenelements have an important role to play in identifying such linearizations. 
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For an unstructured polynomial P G P m (C™ xn ), Higham et al. [131 [H] provide a recipe for 
choosing a linearization by analyzing condition numbers of eigenvalues and backward errors 
of approximate eigenpairs. For structured matrix polynomials, a recipe for choosing a struc- 
tured linearization has been provided in [3] by analyzing structured condition numbers of 
eigenvalues. With a view to identifying optimal and near optimal structured linearizations 
of a structured matrix polynomials, in this section, we analyze the influence of structured 
linearizations on the structured backward errors of approximate eigenelements. It turns out 
that linearizations which minimize structured backward errors also minimize the structured 
condition numbers. Therefore our results are consistent with those in [3]. We thus provide a 
recipe for choosing structured linearizations of a structured matrix polynomial which minimize 
structured backward errors as well as the structured condition numbers. 

For a ready reference, we briefly review some basic results about linearizations of P G 
P TO (C" X ™), for details, see [2TJ1 IT5I [T2l \2~T\ . For our purpose, it is enough to consider the 
vector space Li(P) given by [ST] 

MP) := {L(A) : L(A).(A m _! <g> /„) = v ® P(A), v G C m }, 

where A m _i := [A' m_1 , A m_2 , . . . , 1] T , (g> is the Kronecker product and v is called the right 
ansatz vector for L. Let v := [vi, V2, ■ ■ ■ , v m ] T G C m be a right ansatz vector. Then the 
scalar polynomial p(x; v) := vix m ~ 1 + V2X m ~ 2 + . . . + v m -ix + v m is referred to as the "v- 
polynomial" of the vector v, see [HI [20]. The convention is that p(x; v) is said to have a root 
at oo whenever v\ = 0. Let L(A) = XX + Ye Li(P) be a linearization of P corresponding to 
the right ansatz vector v G C m . Then for x € C™, the following holds 

||L(A)(A m _ l( g)x)|| 2 = |M| 2 ||P(A)x|| 2 , (14) 
\(A m - 1 ®x) T L(\)(A m - 1 ®x)\ = |A£_ lV | |x T P(A)z|, (15) 
|(A m _ 1 ®a;) H L(A)(A m _ 1 ®x)| = |A^_ x t;| \x H P{X)x\. (16) 

Observe from (TT?]) that (X,x) is an eigenelement of P if and only if (A, A m _i eg) x) is an 
eigenelement of L. Consequently, when (A, x) G C x C" is considered as an approximate 
eigenelement of P, it is natural to consider (A, A m _i ® x) G C x C" m as an approximate 
eigenelement of L and vice- versa. We denote the (unstructured) backward error of (A, A m _i<g) 
x) as an approximate eigenelement of L by T}(\, A m _i ® x, L; v) so as to show the dependence 
of the backward error on the ansatz vector v. Similarly, we denote the structured backward 
by rif I (X,A m -- i ®x,L;v), where M G {2,F}. 

Now suppose that (A, x) G C x C" with x H x = 1 is an approximate eigenpair of P. In 
view of (|14| - ()16p . we only need to consider ansatz vectors v having unit norm. We use the 
inequality 

™±1< ^ <1 (17) 

2m ~ HA^xIblKA, l)|| a " 1 1 

which is derived in (Lemma A.l, |13j). 

Theorem 4.1 Let P G P m (C nx ") be regular and L G Li(P) be a linearization of P corre- 
sponding to the normalized right ansatz vector v. Let (X,x) G C x C™ be such that x H x = 1. 
Set A m _i := [A 1 ™^ 1 , . . . , A, 1] T . Then we have 



m + 1 < ??(A, A m _i <8> x,L;v) < 
2m — 7/(A,x, P) 
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Proof: By (JTJ and ([I3|l. we have 
?7(A,A m _i ®a;,L;u) = - 



|L(A)(A m _ 1 ®a:)|| 2 



IM| 2 ||P(A)z|| 



\(A m -i ®a?)|| a ||(A, 1)|| 2 

|(A m _ 1 ® a;) || 2 || (A, 1)|| 2 
|A TO 2 



||(A m _ 1 ®x)|| 2 ||(A, l)|| a 
^(A,x,P) 



r?(A,x,P). 



||A m -i|| 2 ||(A, 1)|| 2 

Hence by (IT71) the desired result follows. ■ 

Theorem 14.11 shows that as far as the backward errors of approximate eigenelements of 
P are concerned, any linearization from Li(P) is as good as any other provided that the 
linearization is associated to a normalized right ansatz vector. In contrast, restricting L in 
DL(P) (see, [20]), it is shown in [4] that the condition number of an eigenvalues A of P is 
increased at least by S(X, v) and at most by V25(\,v), where S(X, v) :— || A m _i|| 2 /|p(x; v)\, 
see also [15] . 

For a structured matrix polynomials, there exists infinitely many structured lineariza- 
tions, see |12[ 1211 |T8] . For the structures we consider in this paper, we consider structured 
linearization from Li(P). For a ready reference, we summarize in Table [5] the condition on 
ansatz vector for a structured linearization, see [THHI]. ^ ne ma trix E in Table [2] is given by 
E = diag{(-ir-\ (-1)™" 2 , . . . , (-1)°}. 



§ 


Structured Linearization 


ansatz vector 


sym 


sksymm 


v g C m 


skew-sym 


skew- sy mm 


v e C m 


T-even 


T-even 


Ew = v 


T-odd 


Ew = — v 


T-odd 


T-even 


Ew = —v 


T-odd 


Ew = v 


Herm 


Herm 


v e m™ 


skew-Herm 


v G iW m 


skew-Herm 


Herm 


v G iR m 


skew-Herm 


v G M m 


T/-even 


H-even 


Eu = v 


T-odd 


Ew = —v 


Tf-odd 


if-even 


Eu = —v 


TZ-odd 


Ev = v 



Table 2: Admissible ansatz vectors for structured linearizations. 



Recall that rj(X,x, P) < ry| f (A, x, P). Similarly, for a structured linearization from Li(P) 
we have r](\, A m _i ® x, L; u) < ?j|^(A, A m _i ® x, L; u), where u is the ansatz vector. With a 
view to understanding the effect of structure preserving linearizations on the backward errors 
of approximate eigenelements of structured matrix polynomials, in this section we compare 
77(A) x, P) and rj%f (A, x, P) with 77^ (A, x, L). 

Corollary 4.2 Let P G § and L G Li(P) be a structured linearization corresponding to the 
normalized ansatz vector v. Then for M G {2,T}, we have 

77| f (A, A m _i g x,L;v) j m+1 
7](X,x,P) ~ V 2m 
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Proof: By Theorem 14. II we have 



rftf(\, A m _i eg) x, L; v) ??(A, A m _i ®x,L;v) jm+1 



7j(A, x, P) 77(A, x, P) V 2m 

Hence the proof. ■ 

4.1 Symmetric and skew-symmetric linearizations 

For a symmetric matrix polynomial P, any ansatz vector v yields a potential symmetric 
linearization. Recall that an ansatz vector v is always assumed to be normalized, that is, 
|| v || 2 = 1. We now show that structure preserving linearizations of symmetric and skew- 
symmetric matrix polynomials have almost no adverse effect on the backward errors of ap- 
proximate eigenelements. 

Theorem 4.3 Let § be the space of symmetric matrix polynomials and P G S. Let L S Lj (P) 
&e a symmetric linearization o/P ™f/i normalized ansatz vector v. Finally, let (A, i)eCx C™ 
&e such that \\x\\2 = 1. X7ien we /iaue 



m + 1 rj| (A,A m _i (g>x,L;?;) r^(A,A m _i eg)x,L;i>) 

2m " r7|(A,x,P) " rj(A,x,P) 



m+1 < 77KA, A m _i g x, L; t>) _ t?(A, A m _i g x, L; t>) < 
2m ~ r?f (A, x, P) ~~ i}(\,x,P) 



Proof: For the Frobenius norm, by Theorem 13. 21 we have 



//■;,(A.A,„_, ,.L:r) _ V 2 " r "" telf |a;Tr|2 ||A,„||> 



V%(\x,P) ^2\\r\\t-\x T r\2 l|A m -i|| 2 || (A, 1)|| 2 ' 

if F (X , A m _i«^,L;^) > 

r}%(X,x,P) ~ V 2m 



where r := — P(A)x. Hence by (IT71) we have — - — ^-g— ^ — - — > • Next, since 



r|| 2 < \/2||r||| - V"- 1 "' |x T r| 2 < V2||r|| 2 , we have 



?7|(A,A m _i Cg)x,L;u) ??|(A,A m _i (8>x,L;u) /- ||A m || 2 

»7|(A,x,P) " r7(A,x,P) " ||A m _i|| 2 ||(A, 1)|| 2 " 



Finally, by Theorem 13.21 we have structured and unstructured backward errors are the 
same for the spectral norm. Hence the desired results follow from Theorem 14. II ■ 

For skew-symmetric linearizations of skew-symmetric matrix polynomials, we have the 
following result. 

Theorem 4.4 Let § be the space of skew- symmetric matrix polynomials and P G §. Let 
L € Li(P) be a skew- symmetric linearization o/P with normalized ansatz vector v. Finally, 
let (A,x) eCxC" be such that \\x\\ 2 = 1. Tften for M e {2,F} we have 



m + 1 < ??^ f (A, A m _i (8) x,L;v) _ ??(A, A m _i <8>x,L;^) < 
2m " ?7m(A,x,P) ~ ??(A,x,P) 

Proof: By Theorem 13.51 we have 

7?| f (A,A m _i ®x,L;v) _ ?7(A,A m _i ®x,L;t;) 
j 7m(A ) x,P) _ r?(A,x,P) 
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Hence desired result follows from Theorem l4.il ■ 

Thus we conclude that for a symmetric/skew-symmetric matrix polynomial a structure 
preserving linearization automatically ensures that the backward errors of approximate eigenele- 
ments are least affected by the conversion the polynomial eigenvalue problem into a generalized 
eigenvalue problem of larger dimension. Moreover, as shown in [4] this choice also ensures 
that the linearization has a mild influence on the structured condition numbers of eigenvalues 
of the polynomial. 



4.2 T-even and T-odd linearizations 

Now we analyze T-even and T-odd linearizations. Note that a T-even (resp., T-odd) polyno- 
mial admits T-even as well as T-odd linearizations which preserve the spectral symmetry of 
the T-even (resp., T-odd) polynomial. 

Theorem 4.5 Let P 6 P m (C nx ™) be a T-even polynomial and (A, x) eCxC" be such that 
\\x\\2 — 1. Let§> e C Li(P) andS C Li(P), respectively, denote the space of T-even and T-odd 
pencils. Finally, let L e 6 § e (resp., L Q £ § 0/ ) be T-even (resp. T-odd) linearization of P 
with normalized ansatz vector v — ~Ev (resp., v — —Hv). Then for M £ {2,F} we have the 
following. 

?7^(A,A m „i (g) x,L e ;v) 



1. //|A| < 1 then J*$± < /MV ' ^^p;' " ^ < V2. 

2 . 7/ |A| > 1 then < ^ x ' L ° ;t,) < V5. 

J 1 1 - V 2m - r?(A,a;,P) 

Proof: First consider the T-even linearization L e . Then by Theorem 13.71 we have 



'//■ ' A. A,„. ; .r.L,:r) ' v ^ ^ ' V ' ' '" ' 

vi\x,P) WallA^HalKA, 1)|| 

where r := — P(\)x. 



(18) 



Now for |A| < 1, we have ||r|| a < J2\\r\\l + ^p^lf 1 ^- !^] 2 < V2\\r\\ 2 . Hence 

by (fT?) we obtain the desired results for the Frobenius norm. 
Again by Theorem 13.71 we have 



?7 2 e (A,A m _i <%>x,L e ;v) 



9 |A| 2 |A T ,v\ 2 . rp io \ ,, . ,, 



V(\X,P) IMIallA^iHalKA, 1)|| 2 



(19) 



Notice that ||r|| 2 <^\\r\\l + |A| 2 \x T r\ 2 < y/1 + |A| 2 ||r|| a for A e C. Hence by (TTT]) 

we obtain the desired result for the spectral norm. 

Next, consider the T-odd linearization L D . Then by Theorem 13. 101 we have 



r//(A,A m _i ®x,L ;u) 



^2||r||I + (^-l)lg^|^rp] ||A m || a 

, (20) 



?7(A,z,P) ||r|| a ||A m _i|| 2 ||(l,A)|| a 
for A ^ 0. Now for |A| > 1, we have 



\rh < Jsllil + (|A|-2 - l)^-^^ < ^HHb. 
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Hence by (I17l) we obtain the desired result for the Frobenius norm. 
Again by Theorem 13.101 we have 



r)2°(\, A m _i ® x,~L a ;v) 



(21) 



»7(A,a;,P) ||r|| 2 ||A m _i|| 2 || (1, A)| 



||A m _! 

we obtain the desired result follows for the spectral norm 



for A 7^ 0. For |A| > 1, we have ||r|| 2 < \/\\r\\% + \\\- 2 j^^ \x T r\ 2 < V2\\r\\ 2 . Hence by 



Remark 4.6 We mention that the bounds in Theorem \4-5\ also hold when P is T-odd with 
the role of T -even and T-odd linearizations are reversed, that is, by interchanging the role of 
L e and L G we obtain the desired bounds. 

Next, comparing n^{X, x, P) with 77§(A, x, L, v) we have the following result. 

Theorem 4.7 Suppose that the assumptions of Theorem \4-.5\ hold. LetS C P m (C nx ") denote 
the set of T -even polynomials. Then we have 

'•TO < 1 : V^< ^^'f-^ * ^ 

2. If |A| > 1 : < ^°(\^rn-i®x K;v) ^ 

J 1 1 ~ V 4m ~ Ti(X,x,P) ~ 

?7 2 (A,A m _i <g) a;, L ; w) 



m+1 1 < '/ 2 k^^-iw^^u, - 

2m A)|| 2 - ry|(A,a;,P) 

Proof: Note that the upper bounds follow from Theorem 14.51 We now derive the lower 
bounds. 

First suppose that |A| < 1. Then it is easy to see that — II e )(A m )|| 2 < ||n e (A m )||2. 
Hence by Theorem 13.71 we have 



1Mb / ||(J-n e )(A TT QiiF V2\\r\\ 2 

A m || 2 \/ l|n e (A m )||2 - ||A m || 2 - 



^(A,z,P) < t] ^ 

On the other hand, by (fTl?|) we have r/| c (A, A m _i ® a;, L e ; u) > — — rrn - ■ Conse- 

||A m _i|| 2 ||(l,A)||2 

quently, by (117|) we have 



77 2 e (A, A ) > ||A m || 2 > 1 J m+1 



4(X,x,P) ~ V5||A m _i|| 2 ||(l.A)||a " 2 

Next suppose that |A| > 1 and consider the T-odd linearization L D . Then it is easy to 
check that ||(/— n e )(A m )|| 2 < ||n e (A m )|| 2 when m is even and the desired result follows 
by similar arguments as above. Now suppose that m is odd. Then it is easy to see that 
|| (J - lie) (Am) HI = |A| 2 ||n e (A m )||i. Hence by Theorem O we have 

^ (XrP)< HI? / . ll(/-n e )(A m )l" VTTWWrh 
%(A ' a; ' Pj - ||A ro || 2 V + ||n e (A m )||l - ||A m || 2 ' 

Further by (f2"T| we have 773° (A, A m _ x ® x,L a ; v) > — J ^ — rr^- Hence by (fTT|) we have 

||A m _i||2|| (1, A)|| 2 



r? 2 °(A,A ro -i ®x,L ',v) \\Knh > 1 /m+1 



4(X,x,P) " l|A m -i||2||(l.A)||2 " V2||(l, A)|| 2 V m 

This completes the proof. ■ 

For T-odd polynomials, we have the following result. 
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Theorem 4.8 Let E> C P m (C nxn ) denote the space ofT-odd matrix polynomials and P € S. 
Let § e C Li(P) and § C Li(P), respectively, denote the space ofT-even and T-odd pencils. 
Finally, leth e G § e (resp., L Q S § J be T -even (resp. T-odd) linearization o/P ratt normalized 
ansatz vector v = "Ev (resp., v = —Yiv). Then for (A, x) G C x C™ wit/i ||x|| 2 = 1, vie have the 
following. 



If | A | < 1 : < & (A ' ~ T " } < 1 



2. If |A| > 1 : =1 ./s±l < < x ^ m ls euerl and 

1 ' " 11(72, A)|| 2 V ™ " r?f(A,x,P) 

4(X,x,P) 



Proof: By Theorem EM we have % s (A,x;P) = _i_./||r||| + , J'"^^ m, 2 1^1 2 - K 

||A m || 2 y IK-/ - ii e j(A m j|| 2 



is easy to see that 



1 + |A| 2 " ||A m ||| 
with equality holds for odd m. Now, by (f2"Tj) we have 



|A| 2 > ||(j-n e )(A m )||| (22) 



|A,„ i|| 2 ||(l, X)\\J Ml + W2 J^^ ] 



r 1 *°(X,A m „ 1 ®x,L ;v) = - „ \ M ^ \/IHII + |A|- 2 l , IA m M2 I^I 2 

' ||(/-n e )(A m )||^ 



Since by ||22J), n J"n Va^ip > |A|~ 2 with equality holds for odd m, we have 



??2°(A, A m _i (g)0J,L o ;ti) < 
?7|(A,x;P) 



(23) 



Also it is easy to check that ||(I- n e )(A m )|| 2 < ||A m || 2 < \/2||n e (A m )|| 2 whenever |A| < 1. 

;j-n e )(A m ; 
||n c (A m )|| 2 



Consequently we have ^{ i n^A^ls ^ — ^his yields 



r] 2 °(\, A m _i (g) x,L ;v) jm+1 



ri(X,x;P) ~ V 6m 

> 1. If 77 

by HU) and ([171) we have 



Next suppose that |A| > 1. If 771 is even then its obvious that nJ^n ^(A ^ 2 )||? — l^l 2 ' Hence 



\M 2 wk^\x T r\ 2 



7^'(A,A m _i(8>a;,L e ;«) < ||A m || 2 V . 

vl(X,x;P) ~ HA^xIblKl, A)|| 2 ^HrHi + IAPI^rl 2 

Further using the fact ll n e( A m)|l 2 < x + | A |2 we have 



ll(i-n e )(A 



m)\\2 



7/ 2 e (A,A m _i ®x,~L e ;v) I jm + 1 



4(X,x;P) ~ ||(V2, A)|| 2 V m 

On the other hand, if m is odd and |A| > 1 then the lower bound follows from the fact that 

l|n e (Am)|| 2 _ ^ ^ This completes the proof. ■ 
\\(I -U e )(A m )\\i l A l 
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The moral of the story is that for computing eigenelements of a T-even matrix polynomial 
P, it is advisable to solve T-even as well as T-odd linearizations of P and then choose a 
computed eigenpair (A, x) from T-even or T-odd linearization according as |A| < 1 or |A| > 1. 
In contrast, when P is T-odd it is advisable to choose (A, x) from T-even linearization only 
when |A| > 1 and the degree of P is even, otherwise choose (A, x) from T-odd linearization of P. 
This choice ensures that the linearizations have almost no adverse effect on the backward error 
of the computed eigenelement (A, x). We arrived at the same conclusion in [4] by analyzing the 
effect of structure preserving linearizations on the structured condition numbers of eigenvalues 
of the polynomial P. 



4.3 Hermitian and //-even linearizations 

First, we consider Hermitian matrix polynomials. Note that a Hermitian matrix polynomial 
admits Hermitian and skew-Hermitian linearizations both preserving the spectral symmetry 
of the Hermitian polynomial. 

Theorem 4.9 Let P e P m (C" x ") be Hermitian and (A, x) eCxC" be such that \\x\\ 2 = 1. 
Let § G {Herm, skew-Herm} and L € § be a linearization o/P with normalized ansatz vector 
v. If A € R then we have 



m+1 J?f-(A, A m _! ® x, L; v) r- lm + l J?f(A, A m _ x ®x,L;v) 

< u r, ^ < V 2 and 4/ — < n 77 — < 1. 



2m ~ ?$ erm (A,a;,P) ~ V 2m ~ ri% erm (\,x,P) 

The same bounds hold when P is skew-Hermitian. 

Proof: First, suppose that S = Herm so that L is a Hermitian linearization of P. For A 6 t, 
by Theorem 13.121 we have 



nr\\l- l £H£\x H r\ 2 



r^(A,A m _i <8)x,L;v) _ \ 2 W^m 1 1 ||A 



m 2 



r^ m (A,x,P) ^IklH-l^rl 2 HA^-iHalKA, 1)|| 

r] S F (A, A m _i g z,L ; u) 
r?^ erm (A,a;,P) 
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where r := — P(A)x. Hence by (TT71) we have — - — ' —, '—^ — > \ h^r^-. Next, since 



||r|| a < J2||r||| - Ig^JJl^rl 2 < y/2 ||r|| 2 , we have 

?7|(A,A m _i 8^L|ij ?y|(A,A m _i (gix,L;i;) /- ||A m || 2 /r 
^ erm (A,z,P) " v(\x,P) ~ V ||A rn _i|| 2 ||(A,l)|| 2 - V ' 

For the spectral norm, by Theorem 13. 121 structured and unstructured backward errors are 
the same when A€l. Hence the desired results follow from Theorem 14. II 

Finally, since the backward errors are the same for Hermitian and skew-Hermitian pencils, 
the above hounds obviously hold for the case when S = skew-Herm. ■ 

This shows that a structured linearization of a Hermitian matrix polynomial does not have 
adverse effect on the backward errors of approximate eigenelements when the approximate 
eigenvalues are real. On the other hand, when the approximate eigenvalues are complex, 
the structured backward errors are not amenable to easy comparisons. Indeed, under the 
assumptions of Theorem l4.9[ when A £ C \ R by Theorem l3.12l a little calculation shows that 

r?|(A, A m _i <g> x, L; v) < j | ||f|| 2 ^ ijf (A, A m _i 



77(A,x,P) "I/ IklH ??(A,x,P) 
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1 


re A 


t ■ 


r = r h := 





imA 







"1 


— imA 


t ■ 


r = r s := 





re A 





re(A^ l _ 1 vx H P(X)x) 
\m{A% l _ 1 vx H P(X)x) 



re(A^ l _ 1 vx H P(X)x) 
\m(A^ l _ 1 vx H P(X)x) 



when S = skew-Herm. 



when S = Herm, and 



Next we consider linearizations of H-even polynomials. Note that an H-even polynomial 
admits H-even as well as -H-odd linearizations and both have the same spectral symmetry as 
that of the polynomial. For purely imaginary eigenvalues of an H-even or ii-odd polynomial, 
we have the following result. 

Theorem 4.10 Let P £ P m (C nxn ) be H-even and (X,x) eCxC" be such that \\x\\ 2 = 1. 
Let S £ {H -even, H -odd} and L £ § be a linearization ofP with normalized ansatz vector v. 
If A £ iM then we have 



m— 1 



X, L: v) 



< %/2 and 



2m ~ ?^- even (A,x,P) V 2m 

The same bounds hold when P is H-odd. 
Proof: The proof is exactly the same as that of Theorem 



1 < r ? |(A,A rn _ 1 



x, L; v) 



v ^- even (X,x,P) 



< 1. 



and follows from Theorem l3.14l 



This shows that a structured linearization of an H-even polynomial has least influence 
on the backward errors of approximate eigenelements when the approximate eigenvalues are 
purely imaginary. On the other hand, when the approximate eigenvalues are not purely 
imaginary, the structured backward errors of approximate eigenelements are not amenable 
to easy comparisons. Indeed, under the assumptions of Theorem I4.10| when A £ C \ iM by 
Theorem I3.14[ we have 



t?|(A, A m _i ®x,~L\v) 
ry(A,.x,P) 




and 



r?i(A,A r . 



1 x, J 



ri(X,x,P) 




where r := — P(A)x and r = r s when § = H-even, and r — rh when § = H-odd. 

The obvious conclusion that we can draw is that real eigenvalues of a Hermitian/skew- 
Hcrmitian matrix polynomial can be computed either by solving a Hermitian or a skew- 
Hermitian linearization. However, for non real eigenvalues it may be a good idea to solve 
Hermitian as well as skew-Hermitian linearizations and choose an eigenpair (A, x) from Her- 
mitian or skew-Hermitian linearization according as r/j < r s or r s < r^. Similar conclusion 
holds for H- even/ H- odd matrix polynomials. These observations are consistent with those 
made in [4] by analyzing the structured condition numbers of eigenvalues. 



5 Structured pseudospectra of structured polynomials 

Let P be a regular polynomial. For A £ C, the backward error of A as an approximate eigen- 
value of P is given by ?7(A, P) := min{?7(A, x, P) : x £ C n and ||x||2 = 1}- Since ?7(A,a;,P) = 
||P(A)x||/||x||2||A m ||2, it follows that for the spectral as well as for the Frobenius norms on 
C nx ™, we have 

<w(P(A)) 



r?(A,P) := 



I A 



m 2 



Similarly, for M £ {2,^} we define structured backward error of an approximate eigenvalue 
A of P £ S by 

?7m(A,P) := min{r/| f (A,x,P) : x £ C" and ||x|| 2 = 1}. 

In this section, we make an attempt to determine ?7|, f (A,P). The backward errors of approx- 
imate eigenvalues and pseudospectra of a polynomial are closely related. For e > 0, the 
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unstructured e-pseudospectrum of P, denoted by cr e (P), is given by (see [SJ[5]) 
cr e (P) = |J WP + AP): APeP m (C" xn )}. 

II|AP|| M < E 

Obviously, we have cr e (P) = {z e C : T](z, P) < e}, assuming, for simplicity, that oo ^ f e (P), 
see [5j [6]. For the sake of simplicity in this section we make an implicit assumption that 
oo ^ cr e (P). Observe that since rj(A,P) is the same for the spectral norm and the Frobenius 
norm, we conclude that a e (P) is the same for the spectral and the Frobenius norms. Similarly, 
when P 6 §, we define the structured e-pseudospectrum of P, denoted by erf (P), by 

o?(P):= (J WP + AP):APe§}. 

|||AP|||M<e 

Then it follows that erf (P) = {z e C : ^ (A, P) < e}. 

Theorem 5.1 Let § € {sym, skew-sym} and P e 8. Then for the spectral norm, we have 
f7§(A,P) = P) and erf (P) = cr e (P). On the other hand, for the Frobenius norm, we 
have ?jf,(A,P) = \/2?7(A,P) and erf (P) = o- e/ ^(P) when § = skew-sym, and n%(X,P) = 
rj(X, P) and erf (P) = er e (P) when § = sym. 

Proof: For the spectral norm, by Theorem 13.21 we have ^f(A,x,P) = n(X,x,P) for all x. 
Consequently, we have ?7|(A,P) = rj(X, P). Hence the result follows. 

For the Frobenius norm, the result follows from Theorem 13.51 when P is skew-symmetric. 
So, suppose that P is symmetric. Then P(A) £ C" xn is symmetric. Consider the Takagi 
factorization P(A) = UYJJ T , where U is unitary and E is a diagonal matrix containing 
singular values of P(A) (appear in descending order). Set cr := £(n, n) and u :—U (:, n). Then 
we have P(A)u ~ o~u. Now define 

T J T 
A A cruu 

1 1 A-m 1 1 2 

and consider the polynomial AP(z) = Y^jLo z^AAj. Then AP is symmetric and P(X)u + 
AP(X)u = 0. Note that for M S {2, F} we have 

»&(A,P) < |||AP||| M = — ^j- = r?(A, P) and hence cr e (P) = erf(P). 

||A m ||2 

This completes the proof. ■ 

When P is symmetric, the above proof shows how to construct a symmetric polynomial 
AP such that A € cr(P + AP) and |||AP|||m = ?7m(A,P). When P is skew-symmetric, using 
Takagi factorization of the complex skew-symmetric matrix P(A), one can construct a skew- 
symmetric polynomial AP such that A £ cr(P+AP) and ||| AP|||m = ?7m(A, P). Indeed, consider 
the Takagi factorization 

P(A) = C/diag(d 1 ,--- ,d m )U T , 
si 

^ q ' 8 i ^ C is nonzero and \sj\ are singular values of P(A). 

Here the blocks dj appear in descending order of magnitude of \sj\. Note that P(X)U = 
C/diag(<ii, • • • , d m ). Let u := U(:, n — 1 : n). Then P{X)u = ud m = ud m u T u. Now define 



where U is unitary, dj := 



A A.- — - ^ '"' 



T 



Mi 

and consider the pencil AP(z) = J2jLo z^AAj. Then AP is skew-symmetric and P(A)w + 
AP(A)u = 0. Further, we have 

% S (A,P) = |||AP||| 2 = a "j; (P { A)) = „(A,P), ^(A,P) = |||AP||| F = V2 a j^ X)) = V2 V (X,P). 
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We denote the unit circle in C by T, that T := {z E C : \z\ = 1}. Then for the T-even or 
T-odd polynomials we have the following result. 



Theorem 5.2 Let S E {T-even, T-odd} and P G § and m is odd. Then for A € T and the 
Frobenms norm we have rf F {\ P) = V2 rj(X, P) and erf (P) nT = <J e/ ^(P) n T. 

Proof: Let AeT. Then by TheoremOand TheoremEIIOl ?7|(A,aj,P) = y/2 ||P(A)x|| 2 /|| A m || 2 
for all x such that ||x||2 = 1. Hence taking minimum over ||a;||2 = 1, we obtain the desired 
results. ■ 

Theorem 5.3 Let § E {Herm, skew-Herm} and P £ S. Then for M E {2,F} and A E R, we 
have rfifiX, P) = r?(A, P) and erf (P) n R = er e (P) n R. 

Proof: Note that P(A) is either Hermitian or skew-Hermitian. Let (/z, u) be an eigenpair of 
the matrix P(A) such that \fi\ — cr m i n (P(A)) and u H u = 1. Then P(A)u = fiu. Define 

_ A 3 n uu H 
3 IIAJH 

and consider the polynomial AP(z) = £?L ^'AA,-. Then AP e § and A E cr(P + AP). 
Further, we have |||AP|||j\/ = — Tr-p— it— — ■ Hence the result follows. ■ 

||A m ||2 

Theorem 5.4 Let § E {H-even, H-odd} and Pe§. Then for M E {2, F} and A S iR, we 
have r/| f (A,P) = r/(A,P) and erf (P) n i R = cr e (P) n i M. 

Proof: Note for A G iR, then the matrix P(A) is again either is Hermitian or skew-Hermitian. 
Hence the result follows from the proof of Theorem 15.31 ■ 

We mention that the above results can be easily extended to the case of general structured 
polynomials where the coefficients matrices are elements of Jordan and/or Lie algebras. 

Finally, we mention that the partial equality crf(L) PI Q = cr e (L) n Q, for an appropriate 
f2 C C, and the minimal perturbations constructed above as well as in section 3 are expected 
to be key tools for solution of certain structured distance problems, see [7]. For illustration, 
consider an H-even polynomial P. We have seen that the eigenvalues of P have Hamiltonian 
spectral symmetry, that is, the spectrum of P is symmetric with respect to the imaginary axis 
in the complex plane and thus appear in the pair (A, — A). Obviously the spectral symmetry 
degenerates if there are purely imaginary eigenvalues. Often in practice the polynomial P 
is obtained by an approximation of the exact problem. Thus it may be the case that even 
though the exact problem has no purely imaginary eigenvalues but due to approximation error 
the polynomial P may have a few purely imaginary eigenvalues. So, the task is to construct 
a minimal perturbation AP such that P + AP is H-even and that the perturbed polynomial 
P + AP has no eigenvalues on the imaginary axis. On the other hand, it may also be the 
case that all eigenvalues of the exact problem are purely imaginary (e.g. gyroscopic system) 
but due to approximation error the polynomial P has a few eigenvalues off the imaginary 
axis. In such a case the task is to construct a minimal perturbation AP such that P + AP is 
H-even and that all eigenvalues of P + AP are purely imaginary. The minimal perturbations 
so constructed and the partial equality of between structured and unstructured pseudospectra 
so established are expected to be key tools in solving these problems and will be investigated 
elsewhere. 

Conclusion. We have derived backward errors of approximate eigenelements of structured 
matrix polynomials. We have constructed minimal structured perturbations that induce the 
approximate eigenelements as exact eigenelements of the perturbed polynomials. The minimal 
perturbations so constructed are expected to be useful in analyzing the evolution of eigenval- 
ues of structured polynomials under structure preserving perturbations. We have analyzed 
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the influence of structure preserving linearizations on the approximate eigenelements of struc- 
tured matrix polynomials. Also, we have provided a recipe for selecting structure preserving 
linearizations so that the linearizations have almost no adverse effect on the approximate 
eigenelements of the structured matrix polynomials. We have briefly analyzed structured 
pseudospectra of structured matrix polynomials and have shown that partial equality be- 
tween structured and unstructured pseudospectra holds for certain structured polynomials. 
These results are expected to be useful for constructing minimal structured perturbations 
that move eigenvalues of the structured polynomials along certain directions which in turn 
are expected to be key tools for solving certain structured distance problems. 
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